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ABSTRACT 

We consider an inhomogeneous model and independently an anisotropic model of primordial 
power spectrum in order to describe the observed hemispherical anisotropy in Cosmic Mi¬ 
crowave Background Radiation. This anisotropy can be parametrized in terms of the dipole 
modulation model of the temperature field. Both the models lead to correlations between 
spherical harmonic coefficients corresponding to multipoles, / and 1 ± 1. We obtain the model 
parameters by making a fit to TT correlations in CMBR data. Using these parameters we pre¬ 
dict the signature of our models for correlations among different multipoles for the case of the 
TE and EE modes. These predictions can be used to test whether the observed hemispherical 
anisotropy can be correctly described in terms of a primordial power spectrum. Eurthermore 
these may also allow us to distinguish between an inhomogeneous and an anisotropic model. 

Key words: hemispherical power asymmetry - inhomogeneous Universe - primordial power 
spectrum - CMBR polarization 


1 INTRODUCTION 


There currently exist several observations which indicate a po¬ 
tential violation of the Cosmological Principle which states that 
the Universe is homogeneous and isotropic on large distance 


trv. teriksen et alj 2004. 20071: Erickcek et al 

I 2 OO 8 I: Hansen et al.l 

20091: Hoftuftetal, 

2009: Pad et al. 2013 

: Ade et al. 2014al: 

Schmidt & Huil 12013 

: lAkrami et alj l2014h 

m the Cosmic Mi- 


crowave Background Radiation (CMBR) which has generated con¬ 
siderable interest in recent years. Parametrization of hemispheri¬ 
cal anisotropy of CMBR t emperature field can be ob t ained by the 
dipol e modulation model jGordon et al.ll200^ l2007l : IPrunet et al.l 
1 20051 : [Bennett et alj201 ih of the following form: 


AT {n) = f{n) ^1 +AA -tij , 


( 1 ) 


where f{n) is an isotropic field, A the preferred direction 
and A am plitude of dipole modulation. Both WMAP and 
PLANCK ilAde et all l2(114al) data confirm this asymmetry 
signal with high signifi c ance. This effect is absent at high-/ 
jPonoghue & Donoghuel IaOOSI : lUanson & Lewis! l2009h . The 
dipole amplitude A and the direction as observed in the WMAP 
five year data are 0.072 ±0.022 and {l,b) = (224°,-22°) ±24° 
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respectively iHoftuft et alj|2009l) . The PLANCK observations are 
consistent with this result. In order to explain the observed large 
scale anisotr o py m a ny theoretical models h ave been proposed 
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2013: 

Mazumdar & Wang 20131: Cai et alj 20131: Liu et al.l 

2013; 

Chang et al.ll201.3l:lMcdonaldll20l4lrTboshll2014h. 


In this paper we are interested in the possibility that the dipole 
modulation might be related to the primordial power spectrum. A 
fascinating possibility is that it might be related to the early inho¬ 
mogeneous or anisotropic phase of the expansion of the universe. 
Such a phase is perfectly consistent with the Big Bang paradigm 
which postulates that the Universe acquires isotropy during infla¬ 
tion. This phenomenon has been explicitly demo nstrated with in 
the framework of the anisotropic Bianchi Models l lWald[l983) . It 
has been found that most of these models become isotropic very 
quickly during inflation. It has been suggested that the anisotropic 
modes, generated during this early p hase may later re-enter the 
horizon iAluril20l^ : lRath et alj2013a) and lead to observed signals 
of anisotropy. These modes will primarily affect the low / multi¬ 
poles of the CMBR and hence are consistent with the observations. 

We consider an inhomogeneous and, independently, an 
anisotropic model of the primordial power spectrum. We assume 
a simple parametrization of the power spectrum in these models. 
We expect that such a form may be related to a fundamental model 
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which might involve an inhomogeneous metric or noncommuta- 
tive space time or some other model of physics at very high en¬ 
ergy scale. Here we use the fact that the observed anisotropy is 
a small effect and hence we only need to keep its contribution at 
leading order. We assume a simple functional form of the k depen¬ 
dence of this contribution and extract the parameters of this func¬ 
tion by making a fit to the temperature data. These parameters are 
then used to predict the anisotropy in polarization data. In particu¬ 
lar we focus on the EE and TE correlations. The power spectrum 
is assumed to be generated by a scalar inflaton field. Since such 
a field does not generate B modes, we don’t consider correlations 
involving these modes. The signature of hemispherical anisotropy 
in CM BR polarization power spectrum has been consider ed ear¬ 
lier in dChang & Wand|2013bl : lNamiodl20l4 Izareill2014l) . How¬ 
ever our emphasis and results are different. In particular we relate 
the observed anisotropy in temperature data to a primordial power 
spectrum and use the resulting power in order to predict the the sig¬ 
nature in CMBR polarization. These results are so far not available 
in literature. 

We next clarify that if we assume homogeneity then it 
is not possible to generate a dipolar power spectrum required 
for an explanation of the hemispherical anisotropy within the 
standard framework of commutative space-time. However there 
are indications that such a power sp ectrum may be obtained 
within a non-commutat iv e space-time {AkoforetalJ l200^j200^ 
Koivisto_&_Mota|_ 201T; iGroeneboom et alT " 2011 : Jain & R^ 


20l4iKotharr et al.l201^ . Here we shall not be concerned with the 


theoretical subtleties in the calculation of this power spectrum and 
its relationship with the ob served CMBR correl ations. These issues 
are partially addressed in iKothari et al.l ( l2015h . However consid¬ 
erable more work is required before this framework can be prop¬ 
erly understood. In any case the inhomogeneous model provides 
a theoretically sound description of the observed hemispherical 
anisotropy. Hence we expect our results for the case of an inho¬ 
mogeneous model to be very reliable. 


2 THEORY 

The CMBR temperature fluctuation and polarizations are fields on 
a sphere, so they can be decomposed in spherical harmonics 

tsT{n)=Y^aT.imVim{n), ( 2 ) 

Im 


E{n) 



(^) 


(3) 


where axjm are the spherical harmonic coefficients and X is the in¬ 
dex denoting temperature (T) or E-Mode polarization (E). In terms 
of primordial density perturbation S(k), they are defined by the 
equation: 

axM = i-i)' To^Tt jd\8{k)Xx (*) (4) 

where To is the present conformal time, Tq is the mean CMBR tem¬ 
perature and Ax {I, k, Tq) is the standard trans fer function employed 
in C AMB. Using the standard notation jZaldarriaga & SeliakI 
Il997h . it can be written as. 


Ax {l,k,ro) 


fj” drS^j'' {k, r) ji (x), 


X = T 
X = E. 


(5) 


The two point correlation function in the presence of dipole 
modulation term in the temperature field can be written as: 

{aTJmOrj'm') — + {^T,lm‘^TJ'm')d m 1 ( 6 ) 

where {aimCi*i„,i)iso = Here Cj is the isotropic angular 

power spectrum and the second term on right is the contribution be¬ 
cause of the dipole modulation t erm in the tempera ture field. It has 
been shown dPrunet et alJl200^ : lRath & Jainiboijh that the dipole 
modulation leads to a correlation between the / and / ± 1 multipoles 
of the temperature field. 

Using Eq.|4l we can write the two point correlations of Eq.[6] 
in terms of the two point correlation function of the primordial den¬ 
sity perturbations {5(k)d*{k')), i.e. the primordial power spectrum 
P(k). In this paper, our primary objective is to consider an inhomo¬ 
geneous primordial power spectrum model and an anisotropic one 
which would give I and / ± 1 correlation explaining the hemispher¬ 
ical anisotropy. In real space, the two point correlation function of 
the density fluctuations 8{x) can be expressed as, 

F{R,X) = {8ix)S{x')), (7) 

where R = x — x' and X = (x + x') /2. If E{R,X) depends only on 
the magnitude R = |R|, then in the Fourier space, the two point 
correlation of 8(k) lead to the standard isotropic power spectrum 



For an anisotropic model E{R,X) depends on both the magnitude 
and direction of R but is independent of X and for an inhomoge¬ 
neous model F{R.X) can’t be independent of X. Next we give de¬ 
tails about the two models. 


2.1 Inhomogeneous Model 

Following ICarroll et alJ 1 I 2 OIOII and iRath et al.l ( |2014|) . we consider 
the real space two point correlation function to have a mild depen¬ 
dence on X and propose the following inhomogeneous model of 


=/i(R) + sin 



5 / 2 W, 


( 8 ) 


where A and 8 are the inhomogeneous parameter and phase re¬ 
spectively. The first term on the right is the standard isotropic and 
homogeneous two point correlator. Since X has the dimensions of 
length, we divide it by the present conformal tim43 Tq in order to 
make it dimensionless. This model assumes an inhomogeneity in 
the direction A, centered around some point whose position is fixed 
by the phase 8. The amplitude of inhomogeneity is parametrized 
by / 2 (R). We expect the absolute magnitude of f 2 {R) to be small. 
Furthermore we expect this term to only show mild oscillations 
and, hence, the absolute magnitude of A should not be large com- 
pa red to unit y . This model differs somewhat from that proposed 
in iRath et aD (1201 4) where the authors expanded the inhomoge¬ 
neous term in powers of X-X and kept only the leading order term. 
Both the models lead to correlations between multipoles / and / -f 1. 
However they differ in the nature of their higher order correlations. 


' One could use some other scale but since the conformal time is related to 
the size of the universe, it is reasonable scale to use. 
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The sinusoidal parametrization is useful since it damps the inhomo- 
geneou s term for large va lues of X. This isn’t the case for the model 
used in iRath et alJ (l2014h . For non-zero values of the parameter 8 
in Eq.[^we obtain correlations between multipoles / and / + 2 also. 
There is currently no evidence for such correlations in data. Hence, 
for simplicity, we set 5 = 0 in our analysis. 

The Fourier transform of the real space two point correlator 
gives 

{s{k)S*{k')) = 

X 

where 

and = (k + K )/2. In this model, the second term of Eq.[^ pro¬ 
duces power asymmetry. We parametrize this function as follows: 

( 10 ) 

where go is the amplitude of the inhomogeneous part, a the spectral 
index and, as defined earlier, Tq is the present conformal time. 


o{k)8^ {k-k’'^ - lj^g(k+) 

% \ 'to 


(9) 


2.3.1 Inhomogeneous Model 


We next compute the contribution to the multipole correlations 
due to the inhomogen eous term. The ca lculation proceeds along 
the lines described in iRath et aD l l2014ll . The relevant quantity, 
Ax^inhomo^ is defined in Eq.[T4l It can be expressed as 


Ax .inhomo (i? 0 


jd^kd^k'Xx {l,k)Ax {l',k') 
^Km{k)yi'm' {k')8ik+) 


X 






where k- = k — k^. We express the spherical harmonics 
as, 




where F"’(cos6) are the Legendre polynomials. After F integra¬ 
tion, we can substitute F = fe ± ^. So up to the first order in the 

inhomogeneous parameter A, we can expand, Ax(/^F) and the 
Legendre polynomials, P”’ (cos 6^/) as. 


2.2 Anisotropic Model 

Eollowine ljain & Ratlil j20l4) and iKothari et all (l2015h . we assume 
that the correlation function depends on the direction of R through 
a dot product with the preferred direction unit vector A. The form 
of the power spectrum in real space is 

=/i(R)+A.R/ 2 (R) ( 11 ) 

where fi (R) is the real isotropic and homogeneous two point corre¬ 
lation function and the second term introduces the anisotropy. Tak¬ 
ing a Eourier transform of the anisotropic term we get 

{S{k)5*ik')) = [pUk) + ik-hik)] 5^ {k-k') , (12) 

In this case we use the same parametrization of g (k) given by Eq. 
[To]which was used for the inhomogeneous model. 

2.3 The Correlation Integrals 

We can calculate the correlation function in Eq.|^in terms of trans¬ 
fer function and the power spectrum by using Eq.|4]for We 

obtain, 

i-i)'-'' I' dhd^k'Ax {l,k)Ax (/',F) 
xY,:„{k)Yi,„,{k')(^5Ck)8{k')) (13) 

This correlation can now be written in the following form, 

Ax,model {id') ( 14 ) 

where the subscript (model) refers to the model, i.e. inhomoge¬ 
neous or anisotropic, being used. As explained earlier, the X sub¬ 
script represents the T and E modes. As has been written above, 
Cx,i is just the isotropic correlation. Here we are primarily con¬ 
cerned with the quantity Ax,„! 0 £/e; dd') which relates I and /', / I'- 

We next compute this correlation for the two models under consid¬ 
eration. 


Ax(/',F) 
ff' (cos 6 ^ 1 ) 


g{k+) 


A Id j\ L ^ a dAx{l',k) 

Ax(l ,A)± — cos0^- - - 

To dk 

A , dP'” (cos dk) 

Pi {cos Ok) ± —Slid dk—!^ -^-. 

kTq acosOf^ 

g(A)±Acose,^^. 


(16) 

(17) 

(18) 


Here 9k is the polar angle of the vector k. The angular integration of 
Ainhomo {h I') implies correlations between I and / ± 1 with m! = m. 
Hence for /' = / -f 1 and m' = m, we obtain. 


Ax,inhomo{ld')-Ti{A7tf (2/-f 1)(2/-f 3) 


where the integrals are defined as follows: 

/i (/,/') = {l + 2)X j\dkAx{l,k)Ax{l’,k)g{k) 

/2 (;,/') = X j\^dkAx(l,k)^{Ax{l',k))g{k) (20) 

h{ld') = ^j^lddkAx{l,k)Ax{l',k)^^ 

where we have used the transfer function given in Eq. Eor the 
inhomogeneous term, we assume the form of g{k) given in Eg. 1101 
Since the parameter A arises as a multiplicative factor in all inte¬ 
grals, it is convenient to absorb it in a redefinition of go- Hence we 
set Ago go- The integrals in Eq.|^are evaluated using a suitably 
modified version of CAMB. 


2.3.2 Anisotropic 

We next compute the contribution to the multipole correlations 
due to the anisotr opic term. We follow the procedure used in 
Ijain & RathI ( l2014h for our analysis. Here the relevant quantity is 
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Ax,aniso which is given hy: 

X iJdk {k) Ax (I,k) Ax [l\k) g(k) (21) 

where 


= 0 


= 8m 


(/ + 1 )^—^ I ;2 _ 

(21+1) (21 + 3) ’’’ (21+1) (21-1) 


Here we’ve used 


|C(e,()))})w(e,<(>) 


(/ + 1)^ — irP- 
(21 + 1) (21+ 3) 


5z',/+l + 


/2 — m2 


(21 + 1) (21+ 3) 


cos Srffl = 




We are interested in only 1,1+1 correlations and by using Eq.|5]we 
obtain, 






— ^mm' ^1 




(4?t)"7S 


{! —m^ 

(21+1) (21 + 3) 


X J k^dkAx{l,k)Ax{l',k)g{k) (22) 


This integral is evaluated by making suitab le changes in the stan¬ 
dard Boltzmann code, CAMB jLewidbOOOh . 


3 DATA ANALYSIS 

Our data analysis procedure closely follows the approach taken in 
iRath et al] ( I20l4) . where the authors have done a detailed analy¬ 
sis of the dipole modula tion term. Her e we use the high resolu- 
tion WAMP-ILC 9 year teennett|[20ll) and PLANCK’s SMICA 
jAde et alj|2014bh data. For masking the galactic portion we use 
KQ85 for ILC and CMB-union mask for SMICA. We use full sky 
maps for computing correlations between different multipoles. The 
full sky maps are generated by filling the masked portions of real 
maps by simulated random isotropic CMBR data. The mask bound¬ 
aries are suitably smoothed by downgrading the high resolution 
maps to a lowe r resolution after applying suitable Gaussian beam, 
as explained in iRath et~^ j20l4) . For the case of the WMAP data 
the contribution due to detector noise is also included while gener¬ 
ating the random map. However contribution due to noise isn’t in¬ 
cluded for the case of SMICA, since the corresponding noise files 
are too bulky. For multipole range 2 ^ 64 the de tector noise 

gives a relatively small contribution dRath et al.ll2014h . Hence the 
error induced by neglecting the noise contribution in the case of 
PLANCK data in the above range is expected to be small. 

The dipole modulation amplitude parameter is estimated from 
the downgraded low resolution maps by computing the statistic 
Sh(L), as defined below. We had seen in section that the dipole 
modulation term in the temperature field, Eq.[T] leads to a non-zero 
correlation between I and I + 1 multipoles. Hence we define a cor¬ 
relation function, C/as, 

r' _ ^(^+1) / „ „* \ 

L/,;+l — 21+1 2^ (23) 

where varies with m and is averaged over the range of 

m for every 1. The correlation is being normalized by the standard 


/(/+!) factor. Since the above correlator shows large fluctuations 
for individual multipoles we define a statistic Sh(L) by summing 
over a range of 21 multipoles in three or six bins, 

Sh(L)= £ Q,,+1 . (24) 

1=1, m,, 


where l^m is t he starting value fo r a bin and Imax = Imin + 20 starting 
from 1 = 2. In iRath et al.l (|2014|) . the authors confined their analy¬ 
sis to three bins, 2^1^ 22, 23 + 43 and 44 ^ ^ 64. The 

downgraded low resolution maps were chosen to have N^ide = 32 
for that analysis. Here we extend this analysis to a large range of 
multipoles with / ^ 127 with the corresponding maps of N^ide = 64. 
For the three bins with 65 I ^ 85, 86 ^ + 106 and 107 127 

we use only the WMAP data since contribution due to noise is ex¬ 
pected to be significant in this range of multipoles. 

We next describe our procedure for estimating the parameter 
A in each of the bins. We first compute the statistic Sh{L) in the 
selected bin from the full sky CMBR map for which the masked 
portions have been filled by randomly generated isotropic data, as 
described above. This estimate of Sh{L) is expected to be biased 
due to the filling of masked portions by random isotropic data. We 
correct this bias by simulations. We generate 1000 full sky sim¬ 
ulated maps which have same characteristics as the re al map, in¬ 
cludin g the dipole modulation. The details are given in iRath et al.l 
1 I2OI4I) . Here we briefly review the main steps in generating sim¬ 
ulated data. We first generate a full sky realization of a random 
isotropic CMBR map. We multiply this map with the dipole mod¬ 
ulation factor, (1 +A/L • n), where the parameters, A and A, are set 
equal to the best fit values obt ained by the hemis pherical analysis 
of the WMAP five year map l lHoftuft et al.ll200^ . We next apply 
the same mask on the simulated map as was used for the case of 
real data. The resulting masked regions are filled with randomly 
generated isotropic CMBR data. This process generates a full sky 
realization of a map which has characteristics similar to that of a 
real map including the dipole modulation. For each value of the 
input parameter A used to generate this map, we generate 1000 re¬ 
alizations and estimate the corresponding statistic Sh{L) by taking 
the average over these maps. The standard deviation of Sh{L) over 
these 1000 maps gives an estimate of the error in the statistic. The 
best fit value of A is obtained by matching the statistic for simu¬ 
lated data with that obtained by real data. The procedure also gives 
an estimate of the error in A. This value of A can be used in order 
to determine the bias corrected estimate of Sh (L) for data. We do 
this by generating a full sky realization of the CMBR map, along 
with the dipole modulation corresponding to the best fit value of A. 
The statistic computed from the resulting map is termed 

The parameters of the theoretical model are computed by fit¬ 
ting the data statistic, The theoretical estimate of the 

statistic, , for a particular set of parameters (ce,go) is ob¬ 
tained my making suitable modifications to the standard Boltz- 
m ann code, CAMB, u sing the cosmological parameters determined 
by lAde et alj (12014bll . We use CAMB to compute the quantities, 
^x,inhomo and ^Xfiniso, which are used to compute the statistic. 
The parameters of the two models are estimated by per- 


^theory 


forming a x'^ minimization. The statistic x'^ is defined as 




cdata 


-S‘ 


,theory 


a,go) 




(25) 


In determining the model parameters, we first set a = 0 and deter¬ 
mine the best fit value of go- Next we determine the best fit values 
of both parameters. 
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Maps 

02^22 

23-5'43 

44^64 

WtLC9 

8.74±5.66 

8.22±3.15 

6.10±2.62 

SMICA 

5.36T5.49 

6.23 ±3.35 

6.60±3.16 


Table 1. The values for data statistic 5^““ x 10^ [mK^) in three bins using 
the WMAP nine year (WILC9) and PLANCK (SMICA) data. The direction 
param eters have been set equal to {/,&) = (224°, —22°) ±24° faoftuft et alj 
l2009l) 




Figure 1. The fit to observed data for three bins (upper graph) and six bins 
(lower graph) for WILC9 map. 


The value for the statistic was calculated independently with 
WMAP 9 year data and the PLANCK data. We use the ILC map 
provided hy WMAP, henceforth called WILC9, and the SMICA 
map provide d by PLANCK. W e use the values of Sff‘‘ that were 
calculated in iRath et al.l ll2014h for the case of three bins. Here we 
also extend this analysis to a larger range of multipoles. These val¬ 
ues are provided in Table[T]for 3 bins and Table|2]for 6 bins. The re¬ 
sulting fits for both the cases are shown in Fig.[T] The direction pa¬ 
rameters have been set to be equal to those o btained by hemispher - 
ical analysis of the WMAP five year data in iHoftuft et al.l (l2009h. 
These match closely with the parameters extracted in Rath et ahl 
( 12014|) by studying correlations between multipoles / and l + l over 
the multipole range 2^1^ 64. The difference in the values of the 
statistic in the first three bins in the Table|2]in comparison to those 
in Table[T]is due to the difference in the resolution of the map used. 
The statistic values in Tables[T]and|2]are extracted by downgrading 
the high resolution map to N^ide = 32 and 64 respectively. In Ta- 
blel^we find that the value of the statistic is significantly lower for 
the bins corresponding to I > 64 in comparison to those of lower / 
bins. Hence we find a clear signal that the dipole modulation effect 
decays beyond / = 64. 


4 RESULTS 

The values of the model parameters obtained by the minimiza- 
tion are given in Tabled These deviat e from the values extracted in 
iRath et al.1 f20li) ;|j' ain & Rati] ( 1201 4t) where only the contribution 
due to the Sachs-Wolfe effect was included. We find that for the 
case of three bins x^ is relatively small. Hence the model provides 
a good fit to data. In fact the x'^ value is too small for the case of 
the SMICA map if we allow a^O. This is probably due to the fact 
that we only have three data points to fit and the error in each is 
relatively large. The fit is not found to be as good for the case of 
six bins. However x^ per degree of freedom is still less than one 
and hence a pure power law model for either the inhomogeneous or 
anisotropic model is acceptable. 

The estimated parameters can be used to predict the signal 
for the case of polarization data. We use the parameters extracted 
for the case of the SMICA map in giving our predictions. In this 
case only data in the three bins, 2^1^ 22, 23 ^ ^ 43 and 

44 ^ Z 64 is used in determining the model parameters. These 
provide a somewhat conservative estimate of the dipole modula¬ 
tion. The predictions obtained by using the WILC9 parameters are 
similar to those obtained by using SMICA. However the amplitude 
in this case is found to be systematically larger by about 15%. Our 
results for TT, EE, TE and ET correlations for the case of the inho¬ 
mogeneous model are shown in Figs.|2and[3l In Fig.l^we plot the 
ratio, defined as, 

ratio = (26) 

Q 

where C; is the standard power and is the correlation be¬ 

tween I and / ± 1 multipoles, defined in Eq.|^ In Fig. we show 
the predicted value of C/./+i for TE, ET and EE modes for the in¬ 
homogeneous model. The corresponding values for the anisotropic 
model are shown in Eigs.|4]and|^ We find that the ratio for the TT 
mode is approximately equal to 0.015 at / = 2 and decays for larger 
values of 1. This is true for both the models. The ratio for the EE 
correlations follow the TT mode closely. The results for the case of 
I ^ 64 are shown in the left panels in these figures whereas those 
for larger / values are shown in the right panels. Since the three bin 
fit may not be reliable over the larger range of I values the results 
for Z > 64 may be treated as qualitative estimates. However the re¬ 
sults for Z ^ 64 are quantitatively reliable and may be compared 
with data in order to determine the validity of the inhomogeneous 
or anisotropic model. 

We give predictions separately for the TE and ET correlations 
since we find that, 

■ ( 27 ) 

The difference between the two arises due to contributions from 
the inhomogeneous or anisotropic terms and can be easily extracted 
from Eq. M This difference is particularly interesting in the neigh¬ 
bourhood of Z ~ 52 where the power C; of the TE mode crosses 
zero. Due to this zero in the denominator, the ratio tends to be 
very large in the neighbourhood of such points. We find that for 
the inhomogeneous model, the TE mode remains positive whereas 
the ET mode becomes negative for Z ~ 40. Hence the ratio for the 
TE mode is positive (negative) for Z < 52 (Z > 52). The ET mode 
shows the opposite trend in the neighbourhood of Z ~ 52. This ap¬ 
pears to be a rather interesting signature of the inhomogeneous 
model which should be testable with the polarization data. Further¬ 
more the anisotropic model (Fig|4ll shows a more pronounced trend 
which is opposite to what is seen for the case of inhomogeneous 
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02^22 

23->43 

44 -> 64 

65 -> 85 

86 ^ 106 

107^ 127 

9.68 ±5.72 

6.78 ±3.39 

6.70±3.13 

1.38±2.91 

4.20 ±2.55 

3.22 ±2.08 


Table 2. The values for statistic 5^“'“ x 10^ {inK^) in 6 bins from WMAP (WILC9). 


No of Bins 

Maps 

Model 


a=^0 


a = 0 





go 

a 


go 


3 

WILC9 

Inhomogeneous 

0.91 ±0.24 

0.38 ±0.08 

0.48 

0.18±0.05 

1.33 

3 

WILC9 

Anisotropic 

I.0I±0.I9 

0.42 ±0.11 

0.50 

0.18±0.05 

1.44 

3 

SMICA 

Inhomogeneous 

0.67±0.2I 

0.32 ±0.07 

6.4 X 10^^ 

0.18±0.06 

0.46 

3 

SMICA 

Anisotropic 

0.75 ±0.23 

0.36 ±0.20 

9.5 X 10-3 

0.17±0.05 

0.48 

6 

WILC9 

Inhomogeneous 

3.25 ±0.89 

0.87 ±0.08 

3.34 

0.04 ±0.02 

7.91 

6 

WILC9 

Anisotropic 

4.13±2.81 

0.94 ±0.24 

3.23 

0.04 ±0.01 

8.05 


Table 3. The model parameters obtained for the case of 3 and 6 bins for the inhomogeneous and anisotropic models. For the case of 3 bins we give results both 
for WILC9 and SMICA. For the case of 6 bins, we give results only for WILC9 for reasons explained in text. 


model. In this case it is the ratio of the ET mode which stays posi¬ 
tive for I < 52 and becomes negative for I > 52 with the TE mode 
showing the opposite trend. Hence we find a qualitatively different 
behaviour for these two models which might be useful in order to 
distinguish between them. 


5 CONCLUSION 

We have investigated the implications of the hemispherical power 
anisotropy within the framework of an inhomogeneous and in¬ 
dependently an anisotropic Universe model. The Universe is as¬ 
sumed to be inhomogeneous and/or anisotropic at very early pre¬ 
inflationary times. This produces a modification of the primordial 
power spectrum which is parametrized in terms of a simple index 
a and amplitude go of ths inhomogeneous or the anisotropic term 
by making a fit to the observed data. The resulting parameters are 
used to make predictions for the polarization data, i.e. EE, TE and 
ET correlations between multipoles corresponding to m' = m and 
I' = l+\. 

In our analysis we mostly emphasize the results for I ^ 64 for 
which the contribution due to detector noise is negligible. For the 
case of WMAP data, however, we also extend our study to larger I 
values including the contribution due to detector noise. We find that 
the signal of dipole modulation is smaller for / > 64 in comparison 
to the lower I values. We use the theoretical parameters extracted 
by fitting data for the three bins corresponding to / 64 for the case 

of SMICA map in showing our predictions. We find that the results 
for the case of WILC9 map are very similar. However the ampli¬ 
tude in this case is found to higher by about 15% in comparison to 
that found in the case of the SMICA map. We separately show our 
predictions for / ^ 64 and for larger I values. For I ^ 64 our pre¬ 
dictions are quantitatively reliable whereas for larger I values they 
should be treated as qualitative estimates. 

We show the results for the ratio of (which is the cor¬ 

relation between / and / + 1 multipoles) to the power C;. We find 
that the results for EE mode follow the TT mode closely. The ratio 
is found to be roughly around 0.02 for 1 = 2 and decays for larger I 
values. We find that the results for the TE mode differ from those of 
ET mode. The difference is emphasized in the inequality shown in 


Eq-|27] For larger I values, in the neighbourhood of Z ~ 52, the two 
modes show qualitatively different behaviour which may provide a 
very clean signature of these primordial models. These predictions 
can be tested in the new data available from PLANCK. Alterna¬ 
tively it may provide stringent constraints on these models. 
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Figure 2. The predicted ratio. Eg. 1261 for the case of the inhomogeneous model for the EE, TE and ET correlations as a function of the multipole /. The TT 
congelations are also shown. The parameter values used for these plots are obtained by fitting the inhomogeneous model to the SMICA map for the case of 
three bins. The left and right plots are for the range / G {2, 3, ... 64} and / G {65,66, ... 300} respectively. The spikes in ET and TE correlations arise because 
the power crosses zero at the corresponding values. 
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Figure 3. The predicted coiTelations, Q^z+i, as a function of the multipole I for the TE, ET and EE modes for the inhomogeneous model. 
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Ratio Of TT and EE Correlations 




Ratio TE and ET Correlations 




Figure 4. The predicted ratio. Eg. 1261 for the case of the anisotropic model for the EE, TE and ET correlations as a function of the multipole /. The TT 
correlations are also shown. The parameter values used for these plots are obtained by fitting the anisotropic model to the SMICA map for the case of three 
bins. 
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Figure 5. The predicted con'elations, C/,/+i, as a function of the multipole I for the TE, ET and EE modes for the anisotropic model. 
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